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A Relativistic 0(a 2 ) Improved Action for Heavy Quarks 
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We extend the Fermilab formalism for heavy quarks to develop an 0(a 2 ) improved relativistic action. We 
discuss our construction of the action, including the identification of redundant operators and the calculation of 
the improvement coefficients. 



1. INTRODUCTION 

A major source of uncertainty in numerical sim- 
ulations of lattice QCD comes from finite lattice 
spacing effects. Since these effects arise at short 
distances, they can be analyzed using an effective 
field theory, as first proposed by Symanzik [Q. 
We may write the effective Lagrangian (LE£) as 



l 



J' 



(1) 



where Sj = dim[Oj]. In this framework lattice 
spacing effects can be systematically removed by 
adding higher dimensional operators to the lat- 
tice action. The aim of this work is to design a 
lattice action such that the coefficients of O(a) 
and 0(a 2 ) terms in Eq. (|l|) can be reduced. 

For heavy quarks with tuq 3> Aqcd, the mass 
introduces an additional short-distance scale into 
the problem, and the Symanzik formalism must 
be modified to separate the short-distance effects 
of both the lattice spacing and the heavy-quark 
mass from the long-distance physics. The Fermi- 
lab formalism [Q] represents one approach to this 
problem. It allows the coefficients Cj in Eq. ([!]) 
to depend on the quark mass. Then it takes 
the Wilson action [||, modified to allow differ- 
ent coefficients for space-like and time-like oper- 
ators. Improved lattice actions are constructed 
by adding higher-dimension operators, again with 
different coefficients for time- and space-like op- 
erators. Ref. |Q] considers operators up to dimen- 
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sion five, to obtain a lattice action for fermions, 
improved to O(a), and valid for quarks with arbi- 
trary mass. In this work we extend the analysis 
of Ref. |^] to include operators of dimension six 
for 0(a 2 ) improvement. 

2. THE ACTION 

At dimension six we have to consider both 
fermion bilinears and four-fermion operators. For 
these proceedings, we concentrate on determining 
the bilinear interactions. The five dimension-six 
bilinears that satisfy axis-interchange symmetry 
are given in Ref. Without axis- interchange 
symmetry, these operators become the seventeen 
operators listed in Table [I]. (The lower dimen- 
sional operators are given in Ref. Q.) 

Following Refs. PP], we use field transfor- 
mations to expose the redundant directions for 
on-shell improvement. Since we are considering 
dimension-six operators in the action we must in- 
clude dimension-two operators in the transforma- 
tions, and we must also consider gauge field trans- 
formations. Writing 



ip — > e J ip, ip — > ?/>e J r , 
we have 



(2) 



J = aei(p + m)+ adij-D + a 2 e 2 {p + m) 2 
- a 2 e F ^a^F^ + a 2 S 2 (j-D) 2 + a 2 5 B iV-B 
+ a 2 S u [ l0 D n-D}, (3) 

and J is the same but with bars over the param- 
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Table 1 

Dimension-six bilinear interactions that could ap- 
pear in the effective Lagrangian. 



No 


Operator 


Parameter 


1 




C2 + £2 


2 


^{ 7 -D,D 2 }V 




3 




^2 + S 2 


4 




82 - S 2 


5 


^{( 7o .D ) 2 ,7-Z>K> 


S u - ~S~u 


6 




S u + S u 


7 


i>~f[D ,E]ip 


e F + e F 


8 


{p[cx-E,j D }^ 


6f — €p 


9 


-ip{'y-D,a-E}'ip 




10 


ip[a-E, j-D]tp 


* 


11 


${iY,-B, l0 D Q }^ 


Sb +~S~b 


12 


tp{i'E-B,f-D}'tp 




13 


tp\iYl-B,j-D]ip 




14 


ip[i'S-B 1 j D ]'ip 


Sb - Sb 


15 


ipjo(D-E - E-D)i> 




16 


ijj~f(DxB + BxD)ip 




17 






eters. 


For the gauge fields 






A + \a 2 t A {D-E - E-D), 


(4) 




+ g 2 a 2 ej{^ t a ip)t a 




A - 


+ A + \a 2 {e A + 8 A ){DxB 4 


BxD) 



- ±a 2 (e A + S E )[Do,E} (5) 
+ g 2 a 2 (ej+8,j)y n t a iP)t a . 

Following the notation of Rcf. §, the e {8) coef- 
ficients label axis-interchange symmetric (asym- 
metric) operators, so that the axis-interchange 
symmetric analysis of Ref. [|J is recovered when 
all 8 coefficients vanish. The transformations of 
Eqs. (||)-((|) generate the operators listed in Ta- 
ble [l] with coefficients which depend on the pa- 
rameters in the transformations. The last op- 
erator of Table [l|, Tp^iDfip, is not generated by 
the transformations, and hence cannot be elim- 
inated. To remove it one must improve the 



nearest-neighbor lattice derivative. 

Not all the operators are independent from 
each other. Indeed, we have used the identities 

2 7o A)7-£>7oA) = (6) 

-j-[D ,E] + {{j D ) 2 ,7-D}, 
2 1 -D l0 D al -D = (7) 
{j-D,a-E} - { 7o A),(7-£>) 2 }, 

2(7-r>) 3 = {~r-D,D 2 } + { 1 -D,iV-B}, (8) 

to remove the operators shown on the left hand 
side of these equations from Table |]. The interac- 
tion ip[a'E,'y-D]i/} (marked "*" in Table |l|) does 
not obey particle-antiparticlc interchange sym- 
metry so it cannot be in the action. Nine in- 
dependent parameter combinations appear in the 
coefficients of the transformed operators. They 
are e 2 _+ £2, ^2_+ 82, 82 - 82, S u + S u , S u - S u , 
Sb + Sb, Sb - S B , £f + £f and e F — If- 

The choice of redundant operators is not 
unique. Considerations, such as calculational 
convenience (or solving the fcrmion doubling 
problem) play a role in the choice. In this work, 
we want to avoid operators which contain higher 
order time derivatives, as they spoil the good 
properties of the transfer matrix of actions with 
Wilson-like time derivatives. There are nine such 
operators, all of which can be eliminated by the 
field transformations. Table [l| shows which pa- 
rameter combination is used to eliminate a given 
operator. In summary, this analysis reduces the 
original seventeen operators to seven, and we can 
write the lattice fermion action as 

S F = S + S B + S E + a 2 ci J ${j-D, D 2 }^ 

+a 2 c 2 J $"/iD? V 4- a 3 c 3 J ^{iH-B^-D}^ 

+a 2 c A J ' $~iq(D-E - E-D)iji (9) 

+a 2 c 5 J 7/7 • (D x B + B x D)tf) 

+a 2 c 6 J^{-y-D,a-E}^/j + a 2 c 7 Jip[i'£-B,'y-D}i/j. 

where So, Sb and Se are given in Ref. M. 
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3. THE COEFFICIENTS 

Conditions on the coefficients of the improve- 
ment operators are obtained by matching on-shell 
quantities in the lattice theory to their continuum 
counterparts. From the lattice fermion propaga- 
tor we derive the dispersion relation 



cosh E = 1 



( M (p)-l) 2 + (/(p)S + 2c 2 Q) : 



Mp) 



(10) 



where (in lattice units) 
fi(p) = 1 + mo + \r s C,p 2 
f(p) = C - 2c 2 - 2c lP 2 



(11) 
(12) 



Si = sinpj, pi = 2sin(j»i/2) and 2Qi = sin2j»i. 
Expanding the energy in powers of p, we obtain 



E = M X 



2M 2 



1 



3 

2=1 



2 ^2 



(p 2 ) 

8M| 



(13) 



where Mi is the rest mass and M2 is the kinetic 
mass. As in Ref. we impose M2 = m q . For a 
relativistic LE£ we also impose Mi = Mi- The 
improvement condition M 2 = M4 yields a relation 
for the coefficient C\ . Rotational invariance (W4 = 
0) gives a condition for C2. These relations were 
already derived in Ref. Q . 

For the remaining operators, C3 through C7, we 
calculate temporal and spatial matrix elements 
with one gluon exchange at tree-level. Using the 
lattice spinors of Ref. Q , we expand the temporal 
matrix element up to and including 0(p 2 ). The 
lattice matrix element is 



(«(^pOlTO£,P)>lat=^(£',0) 
(p' 2 + p 2 - 2p' .p) n ,,,>:!!>', p j 



8M| 



4M|, 



1- 



u(^0), (14) 



where 
1 



8M% 



Cfi 



2mg(2 + TOo) 2 2to (2 + to ) 

C4 



(1 
1 



mo)' 



1 



2c 4 



4M|, 



4M| 



1 + m 



(15) 
(16) 



The right-hand side of Eq. ( |l4| ) must be matched 
to the continuum matrix element, which has Me 



and Me> replaced with m q . One matching condi- 
tion is to set Me 1 — Me, which requires C4 = 0. 
Another matching condition sets Me = M2, 
yielding a condition on c§ and Ce- 



4C 2 



4c £ C 2 



8c 6 



1 



(17) 

mg(2 + m ) 2 to (2 + to ) 1 + m M| 

Eq. (^) is in agreement with a result from Ref. (|] 
obtained in the Hamiltonian formalism. (The ex- 
pression for Mi and the definition of ce can be 
found in Ref. §.) 

The determination of improvement conditions 
on C3, C5 and C7, requires the calculation of the 
spatial matrix element including terms up to 
0(p 3 ). This is currently in progress. 

4. CONCLUSIONS 

We propose a relativistic 0(a 2 ) improved ac- 
tion for heavy quarks. We use the redundant di- 
rections to eliminate all operators with higher or- 
der time derivatives. As a result, our action keeps 
the Wilson time derivative. We have determined 
the mass dependent coefficients of four improve- 
ment operators at tree-level. The determination 
of all remaining coefficients is in progress. 
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